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Motivated by a striking observation of an Fulde-Ferell-Larkin-Ovchinnikov (FFLO) vortex state 
in the heavy fermion material CeCoIns in fields perpendicular to the superconducting planes (H || c), 
superconducting phase diagrams including an FFLO state of uniaxially anisotropic superconductors 
are systematically studied. In the ballistic limit with no quasiparticle's (QP's) relaxation, the high 
field superconducting state in H || c and in the low temperature limit should be not the FFLO 
' state modulating along H, appeared in CeCoIns, but a different vortex state with a modulation 

perpendicular to the field. It is argued that an enhancement near H C 2(0) of the QP relaxation rate, 
presumably originating from a nonsuperconducting quantum critical fluctuation, in this material is 
the origin of the absence of the latter modulated state and of the strange H || c phase diagram in 
which the FFLO state is apparently different from that in H _L c. 
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The recent discovery [l] of a high field superconducting (SC) state in a uniaxially anisotropic heavy fermion 
superconductor CeCoIns in fields parallel to the SC layers (H _L c) has led to renewed interests in the Fulde-Ferrell- 
^ ■ Larkin-Ovchinnikov (FFLO) state 2]. The identification between an FFLO state and the detected high field phase, 
accompanied by a discontinuous iJ C 2-transition 0], is based on an indication of the strong paramagnetic effect in 

S this material [H, 0] and on a derived vortex phase diagram including the discontinuous iz^-transition Q. Based 
on the conventional picture on FFLO states in the vortex-free Pauli limit, however, the presence of an FFLO 
'-^J . state in CeCoIns seems to be an unexpected event in several respects, although this material seems to have a quasi 
two-dimensional (Q2D) electronic structure. First, one needs to clarify why an FFLO state has appeared in the 
material with a weak uniaxial anisotropy although it has not been clearly observed until recently in strongly 
i anisotropic Q2D materials. In particular, the recent observation of an FFLO state in fields perpendicular to the SC 

layers (H || c) 0,1 is the most striking in this sense, because an FFLO state is conventionally expected not to appear 
in this configuration dominated by the orbital pair-breaking. The feature that a flat FFLO transition curve, usually 
expected in the vortex-free Pauli limit, was seen not inHlc but in H || c remains to be explained |7j- Second, an 
observed pressure-induced extension of the FFLO temperature region Q is apparently inconsistent with the fact that 
an FFLO state has appeared in not a material well described by the weak-coupling model but CeCoIns with strong 
electron correlation. 

In this paper, high field phase diagrams including FFLO vortex states of a superconductor with a Q2D electronic 
structure are systematically examined to explain the striking observations in CeCoIns mentioned above on the same 
footing. Phase diagrams are discussed first in the ballistic limit, where the quasiparticle's (QP's) mean free path Zqp 
is infinitely long, and next by assuming a finite Iqp and a slight change of the shape of Fermi surface (FS) . The phase 
diagrams, in both fields parallel and perpendicular to the layers, of the types realized in CeCoIns are obtained only 
when a finite Iqp is assumed in a system with a moderately large Maki parameter cvm- Inclusion of a finite ^qp is 
motivated by two observations: One is the theoretical fact that, in contrast to the conventional ansatz the ground 
state just below H C 2(0) in the ballistic limit is the FFLO state modulating not along H but in the plane perpendicular 
to H @, 

[Til ]. The other is an experimental result suggesting a strong and nonmonotonous field-dependence of Iqp'. 
Recent two transport measurements (l2l . [l3j have shown that Iqp at much lower temperatures than T C (H = 0) is 
anomalously long in low enough fields (i.e., deep in the SC state) and in much higher fields than H C 2(0), implying 
that it is the shortest near H C 2(0). Actually, the pressure induced extension of the FFLO region Q is convincingly 
explained as a consequence of the pressure dependences of T c and Iqp . 

This paper is organized as follows. In sec. II, a theoretical derivation of formulae needed in numerically deriving 
a phase diagram is explained. Numerical results of phase diagrams and thermodynamic quantities are shown and 
discussed in comparison with experimental ones in sec. III. In sec. IV, the obtained results are further discussed. 
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II. THEORETICAL METHOD 



Our analysis starts from the quasi 2D weak-coupling BCS Hamiltonian with the Zecman energy fisH and a d-wave 
pairing interaction which consists of the following three terms 



Ho = dJ2jd 2 r ± (^(r ± ))^ 



(-iV_L + eA) s 
2m P 



(1) 



I d 2 r± Lf(r ± )^ +1 (r ± )+^| 1 (r ± )^(r J 

a, 3 V 
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2 Ej^ B U^) B ^)- (3) 



Here, S CT j(kj_) = 5Z px ApO"^— P-) a J(P+)j P± = P_l ±kj_/2, j is the index numbering the SC layers, pj_ is the 
component of p parallel to the layers, d is the interlayer spacing, and m e is the effective mass of a quasi-particle. 
Further, A p is the normalized orbital part of the pairing-function which, in the case of d 2 .2_ y 2-pairing, is written as 
\J~2(p 2 x — py) in terms of the unit vector p parallel to the layers. Hereafter, the gauge field A will be assumed to 
consist only of a term expressing the external field H, i.e., we work in the type II limit with no spatial variation of flux 
density, because we are interested mainly in the field region near H C 2- Further, afisH = fisH or —fisH is the Zeeman 
energy. In discussing our calculation results, the strength of the paramagnetic effect under a field in the j-direction 
is measured by the Maki parameter a M ,j = V2H^ mb) (0)/H P (0). Here, H P (0) = vrT c /(V2e 7E fi B ) - l.2T c //j, B is the 
Pauli limiting field at T — defined within the weak-coupling BCS model, where 7e = 0.577 is an Euler constant, 
while Hj° lh \o) is the orbital limiting field at T = for fields parallel to the j-direction. 

For the moment, the case with infinite 1 P q (the ballistic limit) will be considered in the weak-coupling approximation, 
and effects of a finite Zqp and strong correlation will be incorporated later. Unless specifically noted, the FFLO state 
on which we focus has a modulation parallel to H. Hereafter, the expressions necessary for exami ning the H-T phase 
diagram in perpendicular fields H || c will be derived by closely following the methods used in Ref.[lll] for H _L c case, 
and hence, just the essential part of the formulation and main results in H || c will be presented below. To describe 
the FFLO state modulating along H || c || z, we can focus on the n = Landau level (LL) modes of the SC order 
parameter A(r) in equilibrium [101 ] . because the n — LL modes are isotropic in nature in the x-y plane and thus, 
cannot accommodate an FFLO modulation in the x-y plane. Then the FFLO state (more precisely, the LO state) is 
described as 

A(r) = y/2T e a e ^ {x, y) cos(Qz), (4) 

where (fo(x,y) is the Abrikosov lattice solution formed in the n — LL under H || z, and Q = Qz is the wavevector 
of the FFLO modulation. Then, the mean field GL free energy density T in n = LL takes the form 

T (n\ 2^(Q) 4 6 
- a a (Q)a e H — a e + — a e 



iV(0)T c 2 ' 2 ~ e ' 3 

= c (0 Ha e ) + c^(a e )q 2 + cS 4 \a e )t (5) 

represented by the amplitude a e and the FFLO order parameter q = 2ttQ^ [5]- Here, N(0) is the density of states 
per spin, = v P / '(2t:T c ) the in-plane coherence length, and Dp the Fermi velocity in 2D case. Microscopic details 
are largely reflected in the expressions of these GL coefficients, a and V m . If necessary, the coefficients ao and V m 
may be expanded in powers of Q 2 : 

oo(Q) = oo(0) + 4V-4V, 

Vi(Q) = V 4 (0)-VPf + V^t, (6) 

where the index "0" of ao indicates the LL index. The ^-dependence of the GL coefficients will be kept up to the 
quartic term so that c^> is given by = al(a^ — a 2 V^ /2). As stressed elsewhere [ll[, inclusion of g-dependence 
of V4 is necessary to keep a stable FFLO state in H 1 c. The same treatment will also be used in H || c. 
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The onset temperature To at which the mean field T/ C 2-transition becomes discontinuous is given as the position at 
which V^{Q m ) becomes negative upon cooling while Vg > 0, where Q m is the equilibrium value of the wave number 
of the FFLO modulation, and a second order transition line HffloCT) is determined as the line on which c^ 2 \a e ) 
becomes negative on cooling, while c^{a c ) > 0. The discontinuous i/^-transition curve is determined by 



3 {V 4 (Q m )) 2 



16 V 6 

Further, by minimizing T with respect both to Q and a e , a 2 is determined by 



OL e {Q m ) 



-Vi{Qm) + ViVijQrn)) 2 ~ 4a (Q m )V 6 
2Vfi 



while 



if a 2 ^ — V^ z> (a e (Q m )) 2 /2 < 0, and q m = otherwise. For instance, in the case with a small but nonvanishing q 2 rll 



q 2 m = (27rQ m Co) 2 



-4 2) + ^ 4 (2) (a e (Q m )) 2 /2 



2(-a 4) + (a e (Q m )) 2 F 4 (4) /2) 



(7) 
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up to 0(q 2 n ), where 
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In numerical calculations we have performed, we always find Sa 2 < 0. That is, the space average of |A| 2 is reduced 
in entering the FFLO state by increasin g H . 

To derive eq.©, the familiar route (El I 111 H3| for deriving a GL action microscopically will be taken. Formally, the 
quadratic term of the GL expression is written as 



where 



K 2 (Tl) 



1 

V 
T 



d 3 rA*(r) 



^ 2 (H) A(r), 



E 



|A p | 2 G £ , ff (p)G_ 



r(-P + n), 



(13) 



(14) 



n 



-id/dr + 2eA(r), and 



ie + oubH — e F 



(15) 



is the quasiparticle Green's function in the normal state in H = 0, where e p is the single particle dispersion measured 
from the Fermi level. Just as in the semiclassical approach 15], the details of FS will be assumed to be reflected just 
in the Fermi velocity vector w and the integral on the FSs when performing the momentum integral. Then, we have 



K 2 = ^V(0)tW|A p 



oo 



isgn(e n ) 



2(ie„ + o^bH) — w • II 
.Y(in / dp/(p)(|A p | 2 exp(iT-Vw-n 
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where 



f(p) 



2nt 



smh(2ntp) \ T ( 



2p B Hp 



(17) 



t = T/T c , and ( )fs denotes the average over FS. 

Similarly, the 4-th order (quartic) term and the 6-th order one of the GL free energy density are written as 



1 

2V 
1 



d 3 rK 4 (U,) A*(n) A*(r 3 ) A(r 2 ) A(r 4 )| rj _ 



T & = — J d J rK 6 (n j )A*(r 1 )A*(r 3 )A*(r 5 )A(r 2 ) A(r 4 ) A(r 6 )| r .^ r , 
where Uj = —id/drj + 2eA(rj). For instance, K4 is given by 

iu = I |Ap| 4 G £ . CT (p)G_ £ ._ ff (-p + nt)G_ e ._ (T (-p + n 2 )G e . (T (p-i-n*-n 2 ) 

= 2nN(0)Ty( - [Sg f £ ' ) \^^ 
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N(0) I n 3 =1 d^ /(]>>) (|A p | 4 exp 
' 3=1 



( piw • n* + /9 2 w • n 2 + p 3 w • n 3 ) 



(19) 



FS 



The corresponding expression of Kq is obtained in the same manner (ll| . 

Applying the parameter integrals used in describing K2 to obtaining the quartic and 6-th order terms [f| [ll[ and 
performing the operation [ll[ exp(ipT c _1 w • II)A(r), we obtain 



a„(0) = -ln(/i)+ / dp 
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where h = H / H ( 2 ™ h) (t=0) = 3.57 e^H, H^"' is the 2D limit of the orbital limiting field Hk° ia ' in H || c, L n (x) is 
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FIG. 1: Cross section in c-a plane of the model FS (solid curves) composed of a corrugated cylinder and a small portion with 
large \w c /w a b\ modelled as a piece of sphere. Any anisotropy in the a-b plane is neglected here. 



the n-th order Laguerre polynomial, 



and 
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r H v F 



(w x +iw y ). 



(21) 



(22) 



To concretely perform the average ( )ps in the above expressions of the GL coefficients, an appropriate FS needs 
to be chosen. To explain the appearance [7| of the FFLO state in H || c, mentioned in Introduction, in CeCoIns with 
a Q2D structure, the use of a simple cylindrical FS is not appropriate: In the case of the purely cylindrical FS with a 
corrugation, the FFLO modulation does not become parallel to H || c [IB] in spite of the field configuration dominated 
by the orbital depairing, reflecting that such a modulation tends to occur in a direction with the largest Fermi velocity 
[17| . This is not in conflict with the presence of the FFLO phase in CeCoIns in H | c because the FS sheet with the 
heaviest mass of quasipaticles in this material is not a pure cylinder with a corrugation but accompanied by a small 
portion with a large \w c /w a b\- See the electron 14-th sheet of FS in Ref.[l8| which has the heaviest effective mass 
and thus, is more effective for a <i-wave superconductivity. In this work, a toy model of FS, sketched in Fig.l, is used 
in which the noncylindrical portion is incorporated as a small piece of the spherical FS with radius hp, the Fermi 
wavenumber in 2D limit. The spanning angle 9 (see Fig.l) measures the size of the noncylindrical portion inducing 
an FFLO modulation parallel to H || c, while the uniaxial anisotropy 7 on of the coherence lengths arises mainly from 
the corrugation of the cylindrical FS. 
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FIG. 2: A typical h v.s. t phase diagram in H || c composed of the mean field transition lines (solid curves) for Iqp = oo, 
«m,c = 6.9, and j an = 2.8. Thin (thick) solid curves imply second order (discontinuous) transition lines for 8 — 2.7 degrees. 
The i/ C 2-curve consists of three solid curves rising upon cooling, while the thin solid curve decreasing upon cooling is -Hfflo(T). 
The open (closed) symbols express the discontinuos i? C 2(T)-transition line (//FFLo(T)-line) for 8 = 4.7 degrees. The I/fflo- 
value for 8 = 4.7 degrees does not decrease unlimitedly but saturates near 0.19 #2o rb) (0) in T-t- limit. The n = 1 LL vortex 
state occurs above the dotted curve. The onset To of discontinuous ii^-transition is indicated by an arrow for each 8. Only 
the -Hfflo(T) curve was sensitive to such a small change of 8- values. 



III. POSSIBLE PHASE DIAGRAMS AND THERMODYNAMIC QUANTITIES 



A typical H || c phase diagram obtained numerically in terms of the tools explained in sec. II is given in Fig. 2, where 
7 a „ is the ratio H^ h \o) / £/i° rb ^ (0) , i.e., the anisotropy of the orbital limited H c2 (0). The figure shows that a drastic 
shrinkage of the FFLO region occurs with decreasing 8 (see Fig.l) reflecting the absence of the H || c FFLO state in 
6 = case [l6j]. The high field ground state just below H C 2(0) in the ballistic limit is formed here not in n — LL but 
in a higher LL, which is n = 1 LL for aM-values of our interest, and has some anisotropic inhomogenuity besides the 
vortices [l^. The n = 1 LL state has a striped structure due to nodal planes || H and can be regarded as another 
FFLO state. In a GL free energy similar to eq.(l) but within the n = 1 LL, its quartic term has a positive coefficient 
near H C 2 for the aM-values of our interest here, and thus, a second order if C 2-transition occurs in t < 0.25 on the 
thin solid line following from a±(0) = (see the first expression of eq. (|20|) ) and rising steeply on cooling. However, 
the presence of the n = 1 LL state is inconsistent with the observations on CeCoIns at low temperatures, in which no 
H C 2 curve rising steeply is seen, and the iJ C 2-transition remains discontinuous even at t ~ 0.015 Q. Note also that 
the change of FS flattens Hyvlo{T) line, while it does not affect the range of the n = 1 LL state and the position of 
the H C 2-Yme. This suggests that other factors than the shapes of FS have to be taken into account to understand the 
H || c phase diagram of CeCoIn 5 . 

It should be stressed that this conclusion does not follow as far as trying to explain only the appearance of the 
FFLO state in H 1 c : As seen in Ref.[13] where the elliptic FS was assumed, the FFLO state in n — LL manages 
to surmount an instability at low temperatures to the in-plane modulated state in n — 1 LL for an appropriate FS 
because the discontinuous H C 2(T) line in H _L c shows a more remarkable rise on cooling than that in H || c. In 
contrast, it is difficult in H || c to, in the ballistic limit, protect the n = LL state from the n = 1 LL one. 

Next, effects of a finite QP's lifetime or Zqp will be considered by neglecting possible h dependences of Iqp for the 
moment. A crucial role of a ii-dependence of Iqp will be pointed out in sec. IV. Possible origins of such a finite Iqp 




FIG. 3: (a) Phase diagram in perpendicular fields (H || c) in the case with a finite mean free path Iqp = 15.5£o- The values 
jan = 2.8, 8 = 4.7 degrees, and «m, c = 6.9 are used. The open circles indicate the curve defined by ai(0) = 0. The arrow 
indicates To. (b) The corresponding one under a parallel field in an antinodal direction following from aM,ai) = 7.69 and the 
same values of Iqp, 8, and y an as in (a). 
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FIG. 4: Field dependences of a dimensionless magnetization m (see the text) in H || c (a) and H 1 c (b) corresponding to 
Fig. 3(a) and (b), respectively. Each arrow indicates the corresponding J/fplo- 



in CeCoIns are an impurity scattering and some magnetic fluctuation accompanying, if any, a quantum critical point. 
Although we focus here on the case of an impurity, the present analysis is qualitatively applicable to the case with 
a magnetic fluctuation as far as the finite QP damping is their main effect on QPs. The primary consequence of 
an impurity scattering is a finite QP's relaxation rate vf/Iqp in anisotropic superconductors satisfying (A p )fs = 0, 
where impurity- induced vertex corrections to a n and V m are negligible 0] . The QP relaxation rate is incorporated in 
eq. (|15[) merely with replacement, 2|e| — > 2|el +«f/£qp, where e is a Fermion Matsubara frequency. In the expressions 
of ea. (|20p . this replacement is represented [f| simply by replacing f(^2 p) there with 

7(Ep)=«p(-^Ep)/(Ep)- ( 23 ) 

A typical example of phase diagrams following from the resulting GL free energy is given in Fig. 3. In H || c, effects 
of the finite Iqp are stronger, and the FFLO state easily shrinks, while the corresponding state in H _L c, as in Fig. 3 
(b), survives over a broad field range keeping a downward transition curve (i.e., dffpFLo/cT > 0). It is a reflection 
of the fact that the H || c FFLO state, supported by a small piece of FS, is fragile and may be easily destroyed by a 
weak perturbation. Such a stronger effect of the finite Iqp on the H || c FFLO state is not surprising once recalling 
that the impurity-induced pinning of vortices in Q2D vortex states is much weaker in the parallel fields. In general, 
as the FFLO state shrinks via a change of FS or a finite Zqp, the Hfflo(T) curve tends to change from a downward 
curve with positive dHp^o/dT to a flat or upward one. More importantly, another modulated state in n = 1 LL 
accompanied by a steep i? C 2-curve was pushed down to T — and was lost by including a finite Iqp in both field 
configurations of Fig. 3. It suggests that a finite QP damping is needed to obtain phase diagrams consistent with those 
of CeCoIns. We note that the ^Qp-value used in Fig. 3 is comparable with the estimated one from thermal conductivity 
data in H > 0.8 (T) P3- 

We note that the above-mentioned disappearance, induced by a finite ^qp, of the n — 1 LL modulated state is a 
consequence of an increase of the LL splitting and thus, cannot be found based on the approach [21] neglecting the 
vortices, i.e., in the limit of vanishing LL splittings. In fact, in contrast to our results [5] and the observation in 
CeCoIns P, 0] , no phase diagrams including a discontinuous _ff C 2-transition are obtained for d-wave pairing systems 
with a finite Iqp in such a limiting case [21( . 

To understand features of the Abrikosov to FFLO transition at -Hfflo, the magnetization M and the specific heat 
in H 1 c have been calculated around -Hfflo- The specific heat jump at -ffFFLO (or T = Tfflo) is given by 

ACfflo _N(0)fdc^\\ AC(0) f dc^ V 

Tfflo (H) 2c( 4 ) \ dt J ~ 20 c( 4 )T c \ dt J ' [ ' 

in terms of the jump value AC(0) at T c in zero field. Calculations leading to Fig. 3(b) show that ACfflo/Tfflo 
is 0.051 AC(0)/T c for t = 0.175 and 0.034 AC(0)/T c for t = 0.075, respectively, which are, up to the factor 2, in 
agreement with the values, (0.065 ~ 0.09)AC(0)/T C , estimated from the data [J, II]. The decrease of ACfflo/Tfflo 
upon cooling (see also Fig. 3 in Ref.[l|) implies that, as Fig. 4 (b) also shows, this transition becomes more continuous 
as the paramagnetic depairing is more effective upon cooling. 

In Fig. 4, the FFLO transitions in the two field configurations are compared with each other through results of the 
normalized magnetization m(T, H) = 8ttk 2 M/(0.12iJ2p J rb ' ) (0)), where k is the GL parameter (i.e, the ratio between the 
penetration depth and the coherence length) defined in low fields. Noting that dT '/da e = dT /dQ — in equilibrium, 



8 



the magnetization M is given simply by 



M = -7V(0)T c 2 a; 



(«o(0))' 

(4 4) )' 



a 2 a 4 



-2 



(25) 



where the prime implies the derivative with respect to H (i.e., (ao(0))' = daa(0)/dH). As mentioned in sec. II, the 
spatial average of |A| 2 decreases due to the appearance of the FFLO modulation. In fact, this decrease of |A| 2 is 
the origin of the slope changes at -Hfflo m thermal conductivity [13] and penetration depth [23[ data. Consistently 
with this | A| 2 -decrease, \m\ also shows an additional reduction, appeared as a kink, on entering the FFLO state from 
below. Note that the kink at -Hfflo in H || c becomes more remarkable rather at lower t in contrast to the tendency 
in H 1 c mentioned above. It implies that a transition to a more fragile FFLO state is sharper, reflecting a rapid 
growth of q near Hfflo, and may become discontinuous [24[ for smaller values of Zqp and/or 9. Thus, for a more 
fragile FFLO state, \m\ decreases more rapidly with increasing H through Hfflo- This feature is consistent with a 
qualitative difference seen between available m-data in H || c and H 1 c [8|,l25||- Further, a sharper change at -/?fflo 
of m tends to reduce the magnetization jump at H C 2- Since a weaker discontinuous iJ C 2-transition should occur closer 
to the virtual second order H C 2 -transition line (i.e., the extrapolations to lower temperatures of the upper thin solid 
curves in Fig. 3), the real _ff C 2-line in such a case will become flatter at lower temperatures. This is a qualitative 
explanation on differences between the H C 2(T) curves in H _L c and H || c. Although it is a conventional wisdom 
that, when the i? C 2-transition is of second order, an FFLO modulation increases the H C 2 value, the coexistence of an 
FFLO state with a flat H C 2 curve in CeCoIns seems to be a consequence of the discontinuous H C 2 transition. 



IV. DISCUSSIONS 



Finally, effects of the electron correlation will be considered here in relation to the p (pressure) dependences of the 
phase diagrams [sj], in which To and Tfflo increasing with increasing p are suggested. The mass enhancement of 
normal QPs, which is a main effect of electron correlation, is incorporated by replacing the Matsubara frequency e 
in ea. lfTB")) by Ze, where Z > 1. Then, by neglecting an e-dependence of Z, it is easily verified that the theoretical 
results in the preceding sections and the figures in sec. Ill, expressed via the normalized field h and temperature t, 
remain unchanged under replacement 

T c = T c (l) -> T C (Z), 
N(0) -> ZN(0), 



Co -> £(Z) 

(2(i! 



ZT C (Z) ' 

A(r) A(r) 



and 



T c ZT C (Z) ' 

a M — > a M {Z) = a M Co t^fS ' > ( 27 ) 

where (ig(Z)H is the Zeeman energy in the case with a mass enhancement included, and the dimcnsionless Fermi 
velocity w/vp is unchanged. The Z dependence in the above A-replacement is familiar in the conventional formulation 
of the strong coupling superconductivity [26[ and implies the difference between the energy gap and the anomalous 
selfenergy of QPs in the SC phase, while the Z dependences of the density of states and the coherence length are 
simply due to the mass enhancement. In general, an enhanced electron correlation should increase hb(Z), (£(Z)) _1 
and thus, om(Z). On the other hand, we have verified that an increase of /xs(Z), as expected, reduces H C 2(0), while 
a decrease of £(Z)results in an increase of H C 2(0). Then, since Z should decrease with p, the p-induced decrease of 
H C 2(0) in H || c p, [27j is mainly a reflection of an increase of £(Z), while the p-induced increase of H C 2(0) in the 
parallel fields 0, [27| is a consequence of p-induced decrease of hb{Z) in H 1 c outweighing the increase of £(Z). 
Then, the p-induced increase of T C (Z) 0] may be one possible origin of the strange p-induced increases 0] of the 
two temperature scales induced by the paramagnetic depairing, Tfflo (H) and the onset To of the discontinuous 
iJ C 2-transition. 

However, the p dependence of the QP damping seems to be a more direct origin of those of the two temperature 
scales if, as suggested in Ref. [1 31 ] , the main origin of the QP damping is a scattering via magnetic fluctuations created 
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by the strong correlation and surviving in t — > limit. In fact, as shown in Fig. 3, we have to assume the presence of 
a small but finite scattering rate ~ ^qp to reach similar phase diagrams to those observed in CeCoIn 5 . According to 
a recent estimation of ^qp based on transport data [15 , [l~3| , Iqp obtained by sweeping h in the low t region relevant 
to the FFLO phenomena seems to be the shortest near H C 2(Q). Then, the results in Figs. 2 and 3 suggest the picture 
that, reflecting a weaker magnetic fluctuation at higher p, the resulting weaker QP damping at higher p would lead 
to increases of the above-mentioned two temperature scales. Further, the strong independence of Iqp(H) suggested 
in Ref . (T3. Jl3l seems to resolve a qualitative disagreement on the -//fflo (T)-line between the experimental phase 
diagrams [ll. foL I22I [23| and Fig. 3(b) in which Iqp was assumed to be iJ-independent : The i?FFLo(^) curve closer 
to the H c2 (T) line in Fig. 3(b) shows a negative curvature (d 2 -Hfflo / dT 2 < 0) in contrast to the experimental one. 
However, if Iqp remarkably decreases with increasing H below H&iT), the part of ^?fflo (T) closer to H C 2(T) should 
be depressed so that -ffFFLo(T) may get a positive curvature. Although discussing a microscopic picture of the 
magnetic fluctuation is beyond the scope of this work, the above arg ument implies that the presence of a magnetic or 
nonsuperconducting quantum critical point near H C 2, suggested [l3l l28j through transport measurements in the high 
field normal state, is the main origin of remarkable differences in the phase diagram from that expected in the weak 
coupling and clean (or ballistic) limit. Then, we speculate that, at higher pressures, the higher LL vortex state with 
a modulation, perpendicular to H, of FFLO type might occur at low enough temperatures in CeCoIns particularly in 
H || c. 

In summary, SC phase diagrams including an FFLO vortex state have been systematically examined to explain 
the H-T phase diagrams of CeCoIn 5 . By examining notable differences in the phase diagrams and thermodynamic 
properties seen between the two configurations, H 1 c and H || c, crucial roles of quasiparticle damping in the FFLO 
state have been pointed out. The origin of differences between conventional theoretical phase diagrams including an 
FFLO state and the observed one of CeCoIns seems to consist in an enhanced quasiparticle damping presumably 
related to a magnetic quantum critical behavior near H C 2- 
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